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We describe the calculation of hydrogenic (one-loop) Bethe logarithms for all states with principal quantum 
numbers n < 200. While, in principle, the calculation of the Bethe logarithm is a rather easy computational 
problem involving only the nonrelativistic (Schrodinger) theory of the hydrogen atom, certain calculational 
difficulties affect highly excited states, and in particular states for which the principal quantum number is much 
larger than the orbital angular momentum quantum number. Two evaluation methods are contrasted. One of 
these is based on the calculation of the principal value of a specific integral over a virtual photon energy. The 
other method relies directly on the spectral representation of the Schrodinger-Coulomb propagator. Selected 
numerical results are presented. The full set of values is available at arXiv.org/quant-ph/0504002. 

PACS numbers: 12.20.Ds, 31.30Jv, 06.20Jr, 31.15.-p 



I. INTRODUCTION 

The evaluation of the Bethe logarithm, in 1947 Jll, was carried out using one of the first automatized devices for the imple- 
mentation of numerical calculations in physics. Today, the evaluation of the (one-loop) Bethe logarithm for the ground state 
of hydrogen to about 10 figures of accuracy can be carried in less than a second on a modern workstation. Consequently, 
one might be tempted ask why there should be yet another paper on Bethe logarithms in the first place? The answer is three- 
fold: (i) In the context of recent efforts toward an improved understanding of the hydrogen and deuterium spectra |01 (see also 
physics.nist.gov/hdel), we have striven to increase the number of states for which this basic quantum electrodynamic correction 
is known, (ii) The Bethe logarithm has been found to follow a characteristic asymptotic structure, expressible in terms of a series 
in inverse powers of the principal quantum numbers. This asymptotic structure has been found to be applicable to wide classes 
of quantum electrodynamic effects in atoms 13 IH. Consequently, it appeared to be of interest to verify these asymptotic 
properties by explicit calculations of the Bethe logarithm for very highly excited states, (iii) The most comprehensive collection 
of Bethe logarithms recorded so far in the literature |0] extends up to the principal quantum number n = 20. Here, we consider 
levels up to n ~ 200. This should be contrasted with recent experimental investigations [7] that were carried out with states of 
principal quantum numbers as high as n = 30. 

Freq uency combs can lead to tremendous simplifications for high-precision spectroscopic experiments (e.g. Ref. IsLl^flol 
fllLfl2l[l^[l4irr^l '). In the future, it should become feasible to carry out high-precision experiments on transitions much more 
effectively than in the past. Rydberg states with long natural lifetimes are rather promising candidates for precision metrology, 
and one might imagine either direct transitions among Rydberg states or optical transitions from, e.g., the metastable 2S state to a 
highly excited D state. Such measurements might contribute to future advances in our knowledge of the hydrogen and deuterium 
spectra. On the theoretical side, the method of least squares |2] allows for a self-consistent adjustment of fundamental constants 
such as the Rydberg constant and the proton charge radius, using experimental input data from more than one transition. The 
investigation reported here is an element of a project (see also Ref. 0) to enlarge the 2002 adjustment of constants lll6ll so as to 
provide optimal predictions for energy levels not contained in this adjustment, consistent with the values of the constants used 
in that adjustment. 

Thus, we here discuss the evaluation of Bethe logarithms using two different methods: 

• an integral representation which is based on analytic calculations using the Sturmian representation of the Schrodinger- 
Coulomb Green function, 

• a spectral representation which relies on known results for the transition matrix elements of discrete-discrete and discrete- 
continuum transitions of hydrogen. 

The two methods are found to be suitable for different ranges of principal and angular momentum quantum numbers. While a 
complete account of all previous work on the Bethe logarithm would result in an excessiv ely l ong list of referen ces, it might 
be instructive and appropriate to recall a few previous investigations on this subject fl7l fl8l fl9l l20l l2ll l22l l23l l24l 12a. l26l 12711 . 
Recently, the Bethe logarithm has been re-evaluated, for selected hydrogenic states, in the context of lower-order terms acting 
as preparatory calculations for higher-order relativistic corrections to the self-energy l28l l29l l30ll . The basic equation defining 
the spectral decomposition has been given in Eq. (2a), and the basic equation for the integral representation has been indicated 
in Eq. (2b) of Ref. 12511 . The evaluation for low-lying states of hydrogen, using the two methods, has previously been discussed 
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in Sees. IIB (integral representation) and and Sees. IIC (spectral representation) of the comprehensive Ref. 12^1 . A similar, 
though less comprehensive, comparison of the two approaches had been made previously in Ref. I22I1 . In the context of the 
spectral representation, we recall that in Table II of Ref. fTvtl and Table I of Ref. fisll one may even find results for the particular 
contributions of the discrete spectrum and of the continuum to the Bethe logarithm of selected low-lying states. In general, 
we found that unexpected numerical difficulties affect the calculation of Bethe logarithms for Rydberg states, especially in 
cases where the difference n — I of the principal quantum number n and the orbital quantum number I is large. Here, we 
attempt to enhance both the range of applica bility of the integral representation as well as the spectral representation, by the 
use of convergence acceleration methods 13 ll 13211 for the evaluation of hypergeometric functions that characterize propagator 
matrix elements (integral representation) and the calculation of infinite sums over discrete virtual intermediate states (spectral 
representation). 

(At least) two further methods are available for the evaluation of Bethe logarithms: one of these is based on the determina- 
tion of approximate eigenfunctions obtained using finite basis sets 1 24], which may be combined with a Neville-Richardson 
extrapolation to yield accurate values, thereby decreasing the required number of functions in the basis set. Basis-set methods 
are also used in calculations of Bethe logarithms in helium (see e.g. Ref. ll33ll34ll V A fourth method relies on a discrete-space 
(lattice) evaluation of the radial component of the Schrodinger-Coulomb propagator l35l l36ll . This method is briefly discussed 
in AppendixlAl 

A somewhat special role is played by circular Rydberg states with n — 1 = / = \m\ @|, whose probability density around the 
atomic nucleus approximately has the shape of a rotationally symmetric, "circular" tire (see Fig. 1 of Ref. |5]), but this shape is 
restricted to the highest possible magnetic angular momentum projection. Because the Bethe logarithm does not depend on to, 
we will refer to all states with n — I = 1 as circular states in this article. Circular states have the highest possible I for given n. 
In the context of the current numerical investigation, it thus appears useful to define a "non-circularity" or "angular-momentum 
defect" £ = n — I > 1. The radial hydrogenic wave functions have the structure of an exponential exp[— r/(nao)], where 
ciq is the Bohr radius, multiplied by a polynomial in r with £ terms. The coefficients of this polynomial have an alternating 
sign pattern. Finite sums whose terms display an alternating sign pattern are notoriously problematic in numerical evaluations, 
because their convergence cannot be accelerated with the methods used for infinite series, and the straightforward summation of 
the terms, using multiprecision arithmetic, is often the only practical route to a reliable numerical evaluation. 

This paper is organized as follows. In Sec.[n] we discuss the integral representation of the Bethe logarithm and its application 
to the calculation of levels with small n — I. In Sec. [HI] we discuss the spectral representation. Some brief conclusions are drawn 
inSec.llVl 



II. INTEGRAL REPRESENTATION 



As is customary for quantum electrodynamic bound-state calculations, we use a system of units in which h = c = eo = 1. 
The Bethe logarithm is a low-energy second-order perturbation which is due to virtual states with one photon mode excited, 
and the atom in a virtual state. This can be seen most clearly by going to the Schrodinger-picture representation of the field 
operators 13711 . A detailed discussion of the transition to the Schrodinger picture for the field operators, together with a basic 
application to bound-state problems, is given in Ref. 13811 . The calculation naturally leads to an integral over the virtual photon 
energy which involves a hydrogenic Green function with an argument z = E n — ui, where u is the energy of the virtual photon. 
In the literature, it is customary to set z = — (Za) 2 m/(2 n 2 1 2 ) (see, e.g., Ref. 0). So t is a variable which parameterizes the 
argument of the Green function in terms of a generalized quantum number "n — > n t." Of course, the variable t has nothing to 
do with temporal evolution. Solving for t, we obtain 



t = 



1 



(1) 



2n 2 to 
(Za) 2 m 



Here, to is the electron mass, Z is the nuclear charge, and a is the fine-structure constant. The photon energy can be expressed 
in terms of t as 



1 - t 2 (Za) 2 m 



(2) 



We denote the reference state by \nlm). The relevant matrix element, which involves the Schrodinger-Coulomb propagator, is 
given by 



1 

H -E + uo{t) 



nil 



(3) 
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The integral representation for the Bethe logarithm In ko(n, I), in terms of P(t), is different for S states (I = 0) in comparison 
to non-S* states (I ^ 0), 



, nW „,, ) = ^ ( P.V.)^d^{^F.( f) + A_8, 



--t'8„i-2\n(n) i,,„ . (4) 

The specification of the principal value is necessary in Eq. (0} because of bound-state poles whose residue gives the one-photon 
spontaneous decay width of an excited atomic state. 

As an example, we consider here the AP state for which the matrix element has the form 

1024 1 7 

P„=4,j=i(t) = - 45 (f_i)8( t + i)8 $ ( 4 ' *) t 75 ~ 1700 * 2 + 9954 1 ~ 2112 4t 6 + 14907 1 8 ] 
2i 2 

+ i—, 7TF7 ttit [15 - 30 1 - 60 i 2 + 150 i 3 + 1547 1 + 15956 1 5 

45 (i - l) 8 (t + l) 8 L 

-154368 1 6 - 142420 i 7 + 1166645 1 8 + 357354 i 9 - 2744516 1 10 - 276066 i 11 + 2046129 i 12 ] . (5) 



The "standard hypergeometric" function which occurs in this expression, is encountered in various previous calculations 1281I2S 

t) = 2 F\ 1, —nt, 1 - nt, 




1 -t x 2 



OO 

fe=0 

The convergence of this series representation near t — is problematic, but it can be accelerated effectively using the combined 
nonlinear-condensation transformation (CNCT) described in Refs. I3lll32il . Using this method, we easily obtain the 40-figure 
result 

lnfc (n = 4,Z= 1) = -0.04195 48945 98085 5486710375 94335 2713418570, (7) 



which is consistent with the 24-figure result given in Eq. (64) of Ref. 13211 and with the 27-figure result given in Table III of 
Ref. jH for this state. 

The number of terms occurring in Eq. is not excessive, but it grows rapidly with the angular momentum defect £ = n — I. 
In addition, considerable numerical cancellation can occur from the typically alternating sign pattern of the polynomial that 
multiplies <I>(n, t). For circular Rydberg states with n — l + 1, the analytic expressions obtained for P(t) are most compact, and 
these states can be well treated using the integral representation, up to very high principal quantum numbers. 

We found that, using the integral representation, numerically satisfactory results can be obtained for states with ( = n — I < 5, 
in the entire range n < 200. However, for ( = n — I > 5, the accuracy obtained using this method was not satisfactory, unless 
an excessively accurate multiprecision arithmetic is used in intermediate steps of the calculation. It is still possible to use the 
integral representation for n — I < 20. In this case, one does not have more than 19 bound-state poles to subtract in forming 
the principal value in Eq. @. As an alternative, one may deform the i-integration contour into the complex plane. However, for 
( = n — I < 20, we found the numerical difficulties to be so severe that a different method of calculation appeared to be called 
for. The more complex structure of the wave function with increasing n — lis illustrated in Fig.^ 

However, before we resort to this different method in Sec. [in] below, we briefly dwell on the application of the integral 
representation to circular states with n = I + 1. It is possible to give a general integral representation for the Bethe logarithm of 
a circular state with I = n — 1, n > 1. This representation is a specialization of and reads 

ln*o(n,n-l) = ~ / dt/(t). (8) 
4 Jo 

Here, 

1 _ j.2 T n f 2 00 

'<«>— Sir -5jr^y+gn*.-.o ■ m 
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FIG. 1: (color online.) The number of terms in the radial wave function grows with the principal quantum number n and with the angular 
momentum defect ( — n — I. In order to illustrate this well-known fact, we here plot the radial probability density of the state with quantum 
numbers n = 40, / = 14 and m = 6, in the plane of constant azimuth <j> = 0. Here, do denotes the Bohr radius ao — h/(amc) = 
0.529 177 2108(18) x 10" 10 m [l^l. The Bethe logarithm for the state under discussion reads In k (n = 40, 1 = 14) = -0.418 087 713 x 10" 4 . 



The term T(k, n, t) is given by 



T(k,n,t) 



-in 



1 fl-t 



+ 



2 4n t 2n (\ + t) 

(2n- l)r(2n) k\ \ 1 + t 

4t 3 T(k + 2n + 2) 
3(1 + f) 4 (-1 - k - n + nt) 



2k 



The polynomial V(k, n, t) reads 



(n - 1) t[P{k, n, t)} 2 T{k + 2n-2) 
3n(l-t) 4 (l-k-n + nt) 



V(k, n,t) = (k + n- l)(2n - 1)(1 - tf - 2tk(k - 1) 



(10) 
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The first term on the right-hand side of dlOl has a singularity at t = (n — l)/n, which corresponds to the decay into the lower- 
lying state with principal quantum number n—1 and orbital angular momentum quantum number I = n — 2. Selected numerical 
values for very highly excited hydrogenic states, obtained using the integral representation, are given in Table|l] 



III. SPECTRAL REPRESENTATION 



The integral representation discussed in SecHTIreflects very closely the physics involved in the original problem, by expressing 
the Bethe logarithm as an integral over the energy of a virtual photon. For computational purposes, a different method can be 
more effective, which relies on available analytic results for transition matrix elements of discrete-discrete transitions l40ll . as 
cited in Eq. (63.2) of Ref. l4lll . and for discrete-continuum transitions [see Eq. (6) of Ref. l42ll or alternatively Ref. ll43ll l. Note 
that in Eq. (6)ofRef. ll42il , the argument of the arccot function in the exponential should be replaced according to n'n/n — > n'/n, 
and that, as pointed out in Ref. itfl . the continuum wave functions used in Ref. l42il are normalized to the energy scale, not to 



TABLE I: Values of Rydberg state Bethe logarithms near n ~ 100, for small angular momentum defect £ = n — I. The 
results displayed here can be obtained using both the integral representation as well as the spectral representation of the Bethe 
logarithm, and both methods were used in order to check the consistency of the results. 



In feo (n, n - 


- C) 


C = i 




C = 2 




C = 3 




C = 4 




n — 


100 


-0.583 308 014 x 10" 


-7 


-0.613 877 681 x 10" 


-7 


-0.645 944 796 x 10" 


-7 


-0.679 594 629 x 10" 


-7 


n = 


101 


-0.566 008 997 x 10" 


■7 


-0.595 371 896 x 10" 


7 


-0.626 158 390 x 10" 


-7 


-0.658 448 703 x 10" 


■7 


n = 


102 


-0.549 387 309 x 10" 


■7 


-0.577 602 405 x 10" 


-7 


-0.607 171 570 x 10" 


-7 


-0.638 170 329 x 10" 


■7 


n = 


103 


-0.533 410 121 x 10" 


-7 


-0.560 532 956 x 10" 


-7 


-0.588 944 368 x 10" 


-7 


-0.618 715 502 x 10" 


-7 


n = 


104 


-0.518 046 496 x 10" 


■7 


-0.544 129 416 x 10" 




-0.571 439 188 x 10" 


-7 


-0.600 042 866 x 10" 


■7 


n = 


105 


-0.503 267 262 x 10" 


■7 


-0.528 359 632 x 10" 


-7 


-0.554 620 643 x 10" 


-7 


-0.582 113 531 x 10" 


■7 


n = 


106 


-0.489 044 896 x 10" 


■7 


-0.513 193 296 x 10" 


-7 


-0.538 455 408 x 10" 


-7 


-0.564 890 904 x 10" 


■7 


n = 


107 


-0.475 353 416 x 10" 


-7 


-0.498 601 821 x 10" 


-7 


-0.522 912 079 x 10" 


-7 


-0.548 340 528 x 10" 


■7 


n = 


108 


-0.462 168 279 x 10" 


■7 


-0.484 558 230 x 10" 


-7 


-0.507 961 045 x 10" 


-7 


-0.532 429 945 x 10" 


■7 


n = 


109 


-0.449 466 295 x 10" 


■7 


-0.471 037 051 x 10" 




-0.493 574 371 x 10" 


-7 


-0.517 128 556 x 10" 


■7 


n = 


110 


-0.437 225 533 x 10" 


-7 


-0.458 014 220 x 10" 


-7 


-0.479 725 687 x 10" 


-7 


-0.502 407 501 x 10" 


-7 



TABLE II: Sample values of Rydberg state Bethe logarithms near n ~ 200, < I < 3. All decimal figures shown are 
significant. The values are consistent with a constant limit for lnfco(n, /) as n — > oo for constant I. In contrast to Table|J 
the values displayed here have been obtained exclusively using the spectral representation. Note that the entries are labeled as 
In ko(n, I) in contrast to the notation In ko(n, n — Q used in Tabled 



In feo ( 


rc,Z) 


I = 


I = 1 




/ = 2 




I = 3 




n = 


190 


0.272 266 958 x 10 1 


-0.490 489 444 x 10" 


i 


-0.993 712 588 x 10" 


-2 


-0.355 864 236 x 10" 


-2 


n = 


191 


0.272 266 942 x 10 1 


-0.490 490 025 x 10" 


i 


-0.993 716 023 x 10" 


-2 


-0.355 866 654 x 10" 


-2 


n = 


192 


0.272 266 927 x 10 1 


-0.490 490 596 x 10" 


i 


-0.993 719 406 x 10" 


-2 


-0.355 869 035 x 10" 


-2 


n = 


193 


0.272 266 911 x 10 1 


-0.490 491 159 x 10" 


i 


-0.993 722 737 x 10" 


-2 


-0.355 871 380 x 10" 


-2 


n = 


194 


0.272 266 896 x 10 1 


-0.490 491 713 x 10" 


i 


-0.993 726 017 x 10" 


-2 


-0.355 873 690 x 10" 


-2 


n = 


195 


0.272 266 881 x 10 1 


-0.490 492 258 x 10" 


i 


-0.993 729 247 x 10" 


-2 


-0.355 875 964 x 10" 


-2 


n = 


196 


0.272 266 867 x 10 1 


-0.490 492 796 x 10" 


i 


-0.993 732 428 x 10" 


-2 


-0.355 878 205 x 10" 


-2 


n = 


197 


0.272 266 852 x 10 1 


-0.490 493 325 x 10" 


i 


-0.993 735 562 x 10" 


-2 


-0.355 880 412 x 10" 


-2 


n = 


198 


0.272 266 838 x 10 1 


-0.490 493 846 x 10" 


i 


-0.993 738 649 x 10" 


-2 


-0.355 882 586 x 10" 


-2 


n = 


199 


0.272 266 824 x 10 1 


-0.490 494 360 x 10" 


i 


-0.993 741 690 x 10" 


-2 


-0.355 884 728 x 10" 


-2 


n = 


200 


0.272 266 810 x 10 1 


-0.490 494 865 x 10" 


i 


-0.993 744 687 x 10" 


-2 


-0.355 886 838 x 10" 


-2 



TABLE III: Values of Rydberg state Bethe logarithms near n ~ 200, 100 < I < 103. This table complements Tab.|l| by 
investigating a range of quantum numbers where both the principal quantum number n as well as the orbital angular momentum 
quantum number I are large. The angular momentum defect £ = n — lis also large for all states in this table. As for the entries 
in TableliTl the values are consistent with a constant limit as n — > oo for given I. 



In fc ( 


n,l) 


I = 100 




I = 101 




I = 102 




I = 103 




n = 


190 


-0.108 830 510 x 10" 


-6 


-0.105 112 540 x 10" 


-6 


-0.101 547 415 x 10" 


-6 


-0.981 275 887 x 10" 


-7 


n = 


191 


-0.109 118 840 x 10" 


-6 


-0.105 395 841 x 10" 


-6 


-0.101 825 818 x 10" 


■6 


-0.984 012 208 x 10" 


-7 


n — 


192 


-0.109 403 831 x 10" 


-6 


-0.105 675 864 x 10" 


-6 


-0.102 101 004 x 10" 


■6 


-0.986 716 920 x 10" 


-7 


n — 


193 


-0.109 685 537 x 10" 


-6 


-0.105 952 663 x 10" 


-6 


-0.102 373 023 x 10" 


■6 


-0.989 390 540 x 10" 


-7 


n — 


194 


-0.109 964 013 x 10" 


-6 


-0.106 226 290 x 10" 


-6 


-0.102 641 927 x 10" 


■6 


-0.992 033 569 x 10" 


-7 


n — 


195 


-0.110 239 309 x 10" 


-6 


-0.106 496 796 x 10" 


-6 


-0.102 907 766 x 10" 


■6 


-0.994 646 499 x 10" 


-7 


n — 


196 


-0.110 511 476 x 10" 


-6 


-0.106 764 231 x 10" 


-6 


-0.103 170 590 x 10" 


■6 


-0.997 229 814 x 10" 


-7 


n — 


197 


-0.110 780 565 x 10" 


-6 


-0.107 028 643 x 10" 


-6 


-0.103 430 446 x 10" 


■6 


-0.999 783 983 x 10" 


-7 


n — 


198 


-0.111 046 625 x 10" 


-6 


-0.107 290 080 x 10" 


-6 


-0.103 687 382 x 10" 


■6 


-0.100 230 947 x 10" 


-6 


n — 


199 


-0.111 309 701 x 10" 


-6 


-0.107 548 591 x 10" 


-6 


-0.103 941 443 x 10" 


■6 


-0.100 480 673 x 10" 


-6 


n — 


200 


-0.111 569 843 x 10" 


-6 


-0.107 804 219 x 10" 


-6 


-0.104 192 674 x 10" 


■6 


-0.100 727 619 x 10" 


■6 
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TABLE IV: Values of Rydberg state Bethe logarithms near n ~ 200, for the highest principal quantum numbers and an- 
gular momenta under investigation in this article. The angular momentum defect is small for the states listed in this Table. 
Consequently, the states listed here can be calculated using both the integral as well as the spectral representation. 

\nk (n,l) 1 = 196 Z = 197 1 = 198 1 = 199 

n = 197 -0.753 369 175 x 10" 8 

n = 198 -0.761 387 888 x 10~ 8 -0.741 963 223 x 10" 8 

n = 199 -0.769 316 490 x 10 -8 -0.749 821 211 x 10 -8 -0.730 786 360 x 10" 8 

n = 200 -0.777 156 469 x 10~ 8 -0.757 591 335 x 10~ 8 -0.738 487 630 x 10~ 8 -0.719 832 864 x 10~ 8 



the momentum scale. The latter fact implies that Eq. (3b) of Ref. 14211 receives a correction according to Eq. (3b) of Ref. |4J 
For the transition matrix elements of selected low-lying states into the continuum, one may alternatively use the formulas given 
in Ref. l45ll . but one should be aware of multiplicative correction factors as pointed out below Eq. (1 1) of Ref. fl7ll . 
The spectral representation is based on the following formula for the Bethe logarithm, 



In ko(n, I) 



2(Za) 4 



i=l 



(H s - En) In 



~2\H S -E n \- 




(Za) 2 m 


m 



nlm 



(12) 



Here, Hs is the Schrodinger Hamiltonian 



P 

2m 



Za 



(13) 



This spectrum of this operator, as is well known, has a discrete part (bound states, E n < 0), and a continuous part (continuum 
states, E > 0). Note that the modulus | Hg — E n | is involved in Eq. (I12> . Otherwise, the argument of the logarithm could become 
negative. The specification of the modulus corresponds to the principal value prescription in Eq. 0. Using the commutator 
relation p l = i m [Hs, r l \, one may easily transform Eq. ill 2> into 



In fcn(n, 



r l (H s - E n y In 



2 \H S — E n 
(Za) 2 m 



2(Za) 4 m 

We assume the wave functions of the continuous spectrum to be normalized according to 

{Elm\E'l'm') = 5(E - E') 5 W 8 mm , . 



(14) 



(15) 



The discrete spectrum is normalized according to (nlm\n'l'm') — 5 nn ' Sw 8 mm ' ■ One may then write down a spectral decom- 
position of d!2i . 



In fco(n, I) 



2(Za) 4 



r l (H s - E n f In 



2\Hs- E„ 



(Za) 2 m 



nlm 



E E E E(^ 



2(Za) 4 m ^ ^ ^ ^ 

y 1 n'=0 l'=l±l m'=-V i=l 



E n r In 



2\E n i — E n \ 



(Za) 2 m 



\(nlm \r l \ n'l'm'^ 



v 1 u l'=l±l m'=-l> i=l L 



2 \E' — E n 
(Za) 2 m 



| (nlm \r l | E'l'm 1 ) | 



= B + C. 



(16) 



Here, B is the bound-spectrum contribution, and C stems from the continuum. Using the dipole selection rules, one immediately 
sees that the angular sums over V collapse to only two nonvanishing terms. We assume the wave functions to have the structure 
(r\nlm) = R n i(r) Yi m (8, <j>) and (r\Elm) = TZ,Ei(r) Yi m (8, <fr), where 9 and are the polar angles. In terms of radial integrals, 
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the quantities B and C read 



B = 



n 3 \ I °° 



(Za) 2 m 



dr r 3 i? nZ (r) i?„'(i_i)(r) 



and 



C = 



/ 4-1 00 

n'=0 



21 



2 LB„/ — E n \ 



{Za) 2 m 



I 



2{Za) A m \2l + l J a 



- J dE' {E 1 - E n f In 



dr r 3 ^(r) i?„'(z+i)(r) 



2 \E' — -En 






(Za) 2 m 




L 



drr 3 R nl (r) K E >(i-i)(r) 



(17a) 



/ + 1 f°° 



f In 



\2\E'-E n [ 




L 


{Za) 2 m 







drr 3 R n i(r) 1l E , {l+l) {r) 



(17b) 



The representations ( I17ai and d!7bt involve only radial integrals; therefore Eq. (63.2) of Ref. 14111 and Eq. (6) of Ref. 14211 can 
be directly applied (with the above mentioned correction in the argument of the arccot function). In the calculation of B, the 
CNCT I3UI32H was used in order to accelerate the convergence of the sum over n' . For the evaluation of C, a simple Gaussian 
integration was found to be appropriate after a suitable variable change that maps the interval E' £ (0, 00) onto the compact 
interval (0,1). 

Selected numerical values for very highly excited Rydberg states with principal quantum numbers n < 200, obtained using 
the spectral representation, can be found in Tables llTllIIII and llVI The values listed in Tableware consistent with the asymptotic 
expansions in AppendixlBl as given in Eqs. iB2\ and flB3l >. and in TablefVl(see also Ref. l46ll ). For circular states, the values given 
in Table II Vl confirm the asymptotic expansion in Eq. JB II . Final numerical calculations were done using the high-performance 
computing facilities of the Max Planck Institute for Nuclear Physics in Heidelberg, and using a cluster of IBM Thinkpad mobile 
workstations l47ll . Advantage has been taken of multiprecision libraries I48ll49ll50ll5lll . 



IV. CONCLUSIONS 



We have presented the evaluation of Bethe logarithms for all the 20100 hydrogenic states with principal quantum number 
n < 200. Two methods have been used: the first method involves an integral representation which reflects the physics of the 
underlying phenomenon in a very direct manner, by expressing the Bethe logarithm in terms of an integral over the virtual photon 
energy (see Sec.[n}. The second method relies on a spectral decomposition of the Bethe logarithm (see Sec.lIIH. In that latter 
representation, the two distinct contributions from virtual discrete states and virtual bound states can be clearly distinguished. 
A third method, which has been used in exploratory work, is briefly described in Appendix [Al Selected numerical data are 
presented in Eq. Q and in the Tables HI— II VI The full set of numerical values is available at 15211 . 

Incidentally, we observe that for all states which simultaneously fulfill 150 < n < 200 and I > 150, the virtual bound states 
give by far the dominant contribution to the Bethe logarithm; indeed, we have \B/C\ > 10 10 for these states [for the definition of 
B and C see Eqs. dl7at and dl7bH . This observation is in sharp contrast to lower-lying states, where \B/C\ is typically smaller 
than one. E.g., the ground state fulfills \B/C\ ~ 0.00483. We conclude that the contribution of virtual bound states as compared 
to virtual continuum states is much more pronounced for Rydberg states as compared to lower-lying states, an observation which 
might appear counter-intuitive at first glance. 
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APPENDIX A: LATTICE SCHRODINGER-COULOMB PROPAGATOR 



The computational implementation of the Schrodinger propagator on a discrete lattice has been discussed in Ref. l36ll . It 
can lead to a computationally cheap evaluation procedure, provided that the required numerical accuracy is not excessive. 
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Numerical difficulties grow for higher excited states, in both one- and two-loop quantum electrodynamic problems, for of a 
number of reasons (more terms and in general a more complex structure of the wave function, more bound-state poles along the 
integration contours, more nodes of the wave function, which translate into numerical cancellations, etc.)- While of course any 
computational method is expandable, we have found it difficult to control the accuracy of the Bethe logarithm, calculated using 
a discrete-lattice representation, for both for high n as well as for high £ = n — I. For high n, the wave function extends over 
many Bohr radii, which necessitates an accurate representation of the Schrodinger-Coulomb propagator over an extended grid, 
which is difficult to achieve with a limited grid size. For high (, the difficulties are enhanced due to the oscillations of the wave 
function which necessitate an even more accurate representation on the grid. 

In exploring this way of calculating the Bethe logarithm, we found it instructive, however, to use routines which lead to 
an explicit diagonalization of the Schrodinger-Coulomb propagator on the grid. The accuracy of the lowest virtual-state en- 
ergy eigenvalues is actually satisfactory, while for higher excited virtual states, the eigenvalues depart rapidly from the exact 
Schrodinger solution E n = —(Za) 2 m/ (2n 2 ). This phenomenon has also been observed in the context of basis-set calculations 
of relativistic effects in atoms which rely on £?-spline techniques, see e.g. Ref. l35ll . On a discrete lattice, one can in principle 
only obtain a discrete spectrum, which for the Schrodinger-Coulomb propagator extends into the positive-energy domain. It 
is then possible to use the eigenvalues and eigenvectors directly in order to evaluate the Bethe logarithm [in the sense of the 
bound-state contribution in Eq. dl7aH . This statement remains true although the higher excited virtual states "energies" on the 
lattice depart very much from the true eigenvalues E n = —(Za) 2 m/ (2n 2 ) obtained for a Hamiltonian acting on L 2 (R 3 ). 

Using the lattice representation, it is easily possible to obtain about 9 decimal figures for the IS Bethe logarithm on a lattice 
which extends to 20 Bohr radii and which is comprised of only 200 nodes. However, we have found it difficult to substantially 
enhance the accuracy, for low-lying states, beyond 20 figures, even if quadruple precision is used in the linear algebra libraries. 
This level of accuracy is of course dwarfed by other available methods, for the concrete problem at hand [see Eq. |7J], and 
therefore the lattice representation has not been pursued any further in the current context of Bethe logarithms for Rydberg 
states. However, we re-emphasize here that the lattice representation can lead to a computationally very efficient evaluation 
of matrix elements of the hydrogenic propagator, a property which has become useful in the calculation of other quantum 
electrodynamic effects for lower-lying states with »i<6@ 15311 . 



TABLE V: Numerical values of the limits lnfco(oo,Z) = lim n _,oo lnfco(n, I) for I = 0, . . . , 10. The evaluation proceeds 
according to methods outlined in Ref. the values communicated here are in agreement with and more accurate than those 
obtained for the range I = 0, . . . , 7 in Ref. Efll . 



1 


In fco(oo, I) 





2.722 654 335 


1 


-0.049 054 544 


2 


-0.009 940 457 


3 


-0.003 560 999 


4 


-0.001 663 771 


5 


-0.000 908 042 


6 


-0.000 548 999 


7 


-0.000 356 923 


8 


-0.000 244 981 


9 


-0.000 175 372 


10 


-0.000 129 830 



APPENDIX B: ASYMPTOTICS DERIVED PREVIOUSLY 

Based on an extrapolation of the numerical data of Tab. I of Ref. Q to higher principal quantum numbers, the following 
asymptotics for the Bethe logarithm of circular Rydberg states had been obtained: 

, , „ „ 0.0248208(6) 0.03814(2) 0.1145(5) 0.166(3) 0.22(2) , „ 
I 3 x In k Q {1 + 1,1) ~ -0.05685281(3)+ ^ + ^— ^ + JT~^ ' (B1) 



Bethe Logarithm for D states 



-0.007 



_ -0.008 
CM 
C 

-0.009 

c 



-0.01 



0.1 0.2 

1/n 

FIG. 2: The Bethe logarithms for / = 2 are in excellent agreement with their asymptotic limit asrn do. We plot here the values In fco (n, 2) 
as a function of n -1 . For n = 190, . . . , 200, the values of lnfco(n, 2) are given in TablellTlfsee Ref. l52ll for a complete list of all relevant 
values in the range n < 200). The limiting value at n^ 1 — 0, which is In fco (00, 2) = —0.994 045 690 x 10~ 2 (see also Tablelvl. has been 
evaluated independently according to Ref. I4H . 




Here, terms of order I k with k > 6 are neglected. For S states, the following asymptotics had been obtained on the basis of the 
numerical data listed in Ref. 0]: 

ln W = 0) * 2.72265434(5) + °- 00000 ° (5) + °- 553 f (5) - °- 5 ", 3(5) + - ^ ■ (B2) 

n n 2 n A n 4 n° 

The corresponding expression for P states reads 

ln fc „M = D * -0.0490545(1) + °- 00000 °( 5 ) + °- 2 ° 53 9 °( 15 ) - °« + _ 3(1) , (B3) 

In an apparently not very widely known paper 14611 , a method has been indicated for the evaluation of the Bethe logarithm in 
the limit of infinite principal quantum number n = 00. This method relies on an asymptotic expansion of the radial integrals 
that enter into Eqs. J17ai and dl7b> . in the limit of an infinite principal quantum number of the reference state. The resulting 
integrals are very slowly convergent, but they permit a completely independent evaluation of In fco(oo, I) which does not rely on 
an extrapolation of data available for lower principal quantum numbers [cf. Eqs. JB 1 i . dB2b and IB3H . Thus, the limiting values 
In fco(cx), I) provide a sensitive independent cross-check of the numerical methods employed in the current calculation. 

In Ref. 14611 . the limits In ko(n = 00, 1) have been evaluated for 1 = 0,. . . , 7 . Here, we generalize the treatment to the range 
I = 0, . . . , 10, thereby confirming the asymptotic limits obtained in Ref. 14611 (see Table IW A comparison to the numerical 
values obtained in the current investigation indicates excellent agreement with the asymptotic values as n — > 00 (see Fig.[2]for 
the case I = 2). A full investigation of the asymptotic structure of Bethe logarithms for large n, including a derivation of the 
subleading terms in the expansion in powers of n _1 , would be very interesting in its own right. Such an investigation would be 
facilitated by the availability of accurate numerical data over wide ranges of n and I, as obtained in the current investigation. 
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